Slow-roll freezing quintessence 
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We examine the evolution of quintessence models with potentials satisfying (V' /V) 2 <C 1 and 
V" /V <C 1, in the case where the initial field velocity is nonzero. We derive an analytic approxi- 
mation for the evolution of the equation of state parameter, w, for the quintessence field. We show 
that such models are characterized by an initial rapid freezing phase, in which the equation of state 
parameter w decreases with time, followed by slow thawing evolution, for which w increases with 
time. These models resemble constant-V" models at early times but diverge at late times. Our 
analytic approximation gives results in excellent agreement with exact numerical evolution. 

PACS numbers: 98.80.Cq ; 95.36.+X 



I. INTRODUCTION 

Cosmological data from a wide range of sources in- 
cluding type la supernovac [lj-3E], the cosmic microwave 
background [3, baryon acoustic oscillations la, B , cluster 
gas fractions 3, Q and gamma ray bursts [9[[lol] seem to 
indicate that at least 70% of the energy density in the 
universe is in the form of an exotic, negative-pressure 
component, called dark energy. 

The dark energy component is usefully parameterized 
by its equation of state parameter, defined as the ratio 
of its pressure to its density: 



w = Pbe/p 



DE- 



(1) 



Observations constrain w to be very close to — 1 . For ex- 
ample, if w is assumed to be constant, then — 1.1 w ^ — 
0.9 0,[l|. 

While a cosmological constant (w = —1) remains con- 
sistent with the observations, a variety of other models 
have been proposed in which w is time varying. A com- 
mon approach is to use a scalar field as the dark energy 
component. The class of models in which the scalar field 
is canonical is dubbed quintessence [l4l - fl7| and has been 
extensively studied. (See Ref. for a recent review). 

A related, yet somewhat different approach is phantom 
dark energy, i.e., a component for which w < —1, as first 
proposed by Caldwell [18j . Such models have well-known 
problems fl9l - [22j (however see [23|, 24| for recent attempts 
to construct a stable model) , but nevertheless have been 
widely studied as potential dark energy candidates. 

These models allow considerable freedom in the choice 
of the potential V ((/>), leading to an infinite set of possi- 
ble models and corresponding behaviors for the evolution 
w as a function of the rcdshift z. It would therefore be 
a considerable simplification if one could find an inter- 
esting subset of such models that converged to a single 



trajectory or a well-defined set of trajectories for w(z). 
One approach that yields such a simplification was con- 
sidered by Scherrer and Sen [25} , who examined the evo- 
lution of a scalar field in a "nearly flat" potential, where 
the flatness condition consisted of the slow-roll conditions 
familiar from inflation: 
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As in Ref. [25j, we will assume for definiteness that 
dV/d(j) < and A > 0, but of course none our results 
will depend on this choice. Note that while equations ([2]) 
and ([3]) are the familiar slow-roll conditions from infla- 
tion, the evolution of the scalar field is very different for 
the case of quintessence, since in that case one must also 
include the effect of the matter density on the expansion 
rate. 

Ref. [15 considered models satisfying equations ([2]) 
and ([3]) in which the field is initially at rest, so that 
w ~ — 1 at early times, and w increases at late times 
as the field rolls down the potential; in the terminology 
of Ref. (26[, these are "thawing" models. Then equation 
([2]) ensures that w remains close to —1, while equations 
© and © taken together indicate that A is nearly con- 
stant. For all potentials satisfying these conditions, it 
can be shown that the behavior of w can be accurately 
described by a unique function of 0^ (the fraction of the 
total density contributed by the quintessence field, where 
we assume a flat universe) and the (assumed constant) 
value of A (25[. In 27 1 this result was extended to phan- 
tom models satisfying Eqs. ([2]|3j. 

In this paper, we relax the initial condition that <f> = 
d<p/dt — 0, so that w 7^ —1 initially, but we retain both 
slow roll conditions on the potential (equations [5] and [3]) . 
This allows us to extend the formalism of Ref. [!i| to 
models in which w decreases toward —1 at late times, 
rather than increasing away from —1; the former models 
are called "freezing" models (26[. 

Note that the slow roll conditions, Eqs. (HUH]), while 



2 



sufficient to ensure w ~ — 1 today, are not necessary, 
and many other attempts to ciassify or simplify the set 
of quintessence trajectories have been proposed. For ex- 
ample, if equation @ holds, but equation ([3]) is relaxed, 
one still has w m — 1, but there is now an extra degree of 
freedom, the value of V" /V. Instead of a single solution 
for the evolution of w, one obtains a well-defined family 
of solutions [28| . (This family of solutions includes the 
slow- roll solution of Ref. [25| as a special case in the limit 
where V" /V — > 0). These models were explored in more 
detail in Refs. HaHJ. 

Other attempts to systematize the behavior of dark- 
energy quintessence fields based on other assumptions 
have been given in Refs. [3214391 ]. Where these other 
approaches overlap those of this paper, further discussion 
and comparison with our results will be given below. 

In the next section, we examine the evolution of w 
for slow-roll potentials with the assumption that 4> 
initially. In Sec. Ill, we compare our predictions for 
w(a) with exact numerical results. Our conclusions are 
discussed in Sec. IV. 



(See Ref. [25| for the details of this derivation). 

Consider a scalar field moving in a potential that satis- 
fies equations @ and ([3]) , and assume we are in a regime 
such that 1 + w <C 1. Note that in relaxing the initial 
condition that = 0, it is possible to have transient evo- 
lution with w far from —1 even when equations @ and 
§3§ are satisfied. Here we assume that the evolution has 
proceeded far enough that l + w <C 1, but w ^ — 1, when 
we begin to examine the evolution. With these assump- 
tions, equation © simplifies to [25| 



dw 



2(1 + w) 2 y/3(l + W ) 
'fy(l-fy) 3 "(l-fyWTV 



(9) 



where Ao is the (assumed constant) value of A. The gen- 
eral solution of equation ([9]) is 



l + w= -X 2 



n, 



- i 



tanh 



(10) 



II. EVOLUTION OF w: ANALYTIC RESULTS 

We will assume that the dark energy is provided by a 
minimally-coupled scalar field, </>, with equation of mo- 
tion given by 



3H<j>+^- = 0, 



where H is the Hubble parameter, given by 



(4) 



The constant C parametrizes different initial values of 
(f>, and therefore characterizes different evolutionary tra- 
jectories. Ref. [25[ considered only the case C = 0, which 
corresponds to 4> = initially, giving 



1 



1 



1 I tanh 1 \/TL 



(11) 



H = 



(5) 



Here a is the scale factor, p is the total density, and 
we take SttG = 1 throughout. Equation Q indicates 
that the field rolls downhill in the potential V {(f)), but its 
motion is damped by a term proportional to H. 

The pressure and density of the scalar field are given 

by 



(6) 



and 



Equation (fTTj) . derived as one of the main results in Ref. 
[251] . displays a purely thawing behavior, as the field be- 
gins initially with w = — 1, and w increases as <fi rolls 
down the potential. Note that at early times, when 
fi^ <C 1, we can expand equation (jlip to give 
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V/2 + Aft3/2 1^5/2 

3 15 35 



(12) 



Retaining only the first term in this expansion gives 



i + w = —x n 4 



(13) 



P=\ + V{<t>), 



respectively, and the equation of state parameter, w, for 
the quintessence field is given by equation ([1}. 

As in Ref. [25[, we note that the equations simplify if 
we express w as a function of fi^, giving 



dw 

dflrh 



-3(1 + w)(l -w) + A(l - w)y/3(l + w)n^ 

3(-«;)n^(i-n^) 



This a gree s with the result of Cahn, dc Putter, and 

(7) Linder [35(, who defined a "flow parameter" F = (1 + 
u;)/A§r^ and showed that F — 4/27 in the limit consid- 
ered here. 

Equation (fTU|) when C / corresponds to the case 
where the field has cb ^ initially. The case C < corre- 
sponds to a field rolling down the potential, while C > 
gives solutions for the field rolling uphill initially. While 
the latter possibility might seem somewhat contrived, it 
has been considered previously pol - l42l ] and we include it 

(8) here for completeness. 



3 



The actual value of C depends in the initial value of 
4>, but it is easiest to evaluate C in terms of the value 
of w corresponding to some initial value of 0^. If we 
designate these values as Wi and fi^j, respectively, and 
take an early enough epoch that Q,^ <C 1, then we have 

c = ±y /3(i + m)n^ (w) 
Ao 

where the sign of C is determined by whether the field is 
rolling initially up or down the potential. In terms of 0, 
we simply have 1 + lu, « <j> 2 /Vj, where V* is the (initial 
but assumed nearly constant) value of the potential. 

In equation (| 1 Of) w enters only in the combination (1 + 
if) /Aq. In the left panel of Fig.[TJ we display this quantity 
as a function of £1^ for the indicated values of C < 0. 
Except for C — 0, the generic behavior of these models 
is an initial freezing period, during which w decreases 
rapidly with time, followed by a period of thawing, with 
a much slower increase of w with time. 



III. COMPARISON WITH EXACT EVOLUTION 

We now test the accuracy of the analytic approxima- 
tion derived in the previous section by comparing our 
analytic predictions for w(a) with an exact numerical in- 
tegration of the equations of motion. Although we expect 
our approximation to apply to any potential with a re- 
gion in <f> satisfying equations ^) and ([3]) , we will consider 
here two representive examples: 

• Linear potential 

V (0) = V Q - a<t> (18) 
where take a = 0.1. 

• Power-Law potential 

V{4>) = V cj>- n (19) 
where we take n = 0.1. 



The value of f2^ for which 1 + w reaches a mini- 
mum, corresponding to the transition between freezing 
and thawing behavior, can easily be derived from equa- 
tion fTU|): 



C 



'f2fli — tanh 



(15) 



In both of these cases the parameters of the potential 
were chosen so as to satisfy the slow-roll conditions (equa- 
tions [2] and [3]), so that A is roughly constant and <C 1. We 
set the initial value of w to u>, = —0.95 and Ao = —0.08. 
Then the initial value of fl^, is chosen to give C = ±0.1, 
C = ±0.5 and C = ±1. (Numerical tests of the C = 
case can be found in Ref. (HI)- 



Thus, models with — C £ 0.5 are still freezing at the 
present. On the other hand, it is clear from Fig. 1 that 
models with — C ^ 0.1 are nearly indistinguishable from 
the C = case, as they enter the thawing period when 
^ < 1. 

In the right panel of Fig. [TJ we show the corresponding 
evolution of (1 + u>)/A 2 as a function of il^ for C > 0. 
These models show more complex behavior, as the field 
first rolls uphill, stops (giving w = —1), and then rolls 
back down again. The value of at which w = — 1 is 
given by 



C 



(1-fi) 



tanh 



(16) 



Since we always assume 1 + w <C 1, we can express 
w as a function of the scale factor a by taking Q^, to be 
well-approximated by its value for a ACDM universe [25| : 



1 



%o 1 



-3 



(17) 



Equations (|TU)) and (fT7j) together then provide an ana- 
lytic approximation for w(a) for these models. In Fig. [2] 
we show (1 + w)/Xq as a function of the scale factor for 
the same set of models as in Fig. [I] where we have taken 



n 



00 



0.7 and a = 1 at the present. 



Our results are displayed in Figs. [3]and|U Note that 
the agreement between our analytic approximation and 
the true (numerically-calculated) evolution is excellent. 
This demonstrates that in the regime where the slow roll 
conditions hold, the evolution of w converges to a com- 
mon set of tracks irrespective of the details of the under- 
lying quintessence potential. 

In the limit where A = 0, these models correspond to 
a field evolving in a flat potential, V = Vq; these models 
have been dubbed "skating" models |43j . Since we are 
taking A 1, it is reasonable to ask whether the freezing 
portion of these trajectories can be well-approximated by 
a pure skating trajectory. 

An approximate skating trajectory will be realized 
whenever the first term in the numerator of equation ([8]) 
dominates the second term; this corresponds to the con- 
dition 



A 2 

1 + w > — fl 

O 



<$>■ 



(20) 



Combining this result with equation (flU)) (for C < 0) 
gives the following condition for the evolution to approx- 
imate a skating solution: 



C > tanh' 



(21) 
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FIG. 1: The value of (1 + w)/\q as a function of fi^, for the indicated values of C. Here w is the equation of state parameter 
for the quintessence field, fi^ is the fraction of the total density contributed by the quintessence field, and Ao is the value of 
— (l/V)(dV/d(j>) (taken to be constant). 




FIG. 2: As Fig. [T] with the scale factor a (taken to be 1 at the present) on the horizontal axis, and fi^o = 0.7. 




FIG. 3: Numerical versus analytic evolution of the equation of state parameter, w, for the linear potential. From left to right 
the plots correspond to \C\ = 0.1, \C\ = 0.5 and |C| = 1 respectively. The sign of C is indicated on the figure. Solid curves 
correspond to the exact (numerically-derived) evolution, while the dashed curves indicate the analytical approximation from 
equations (|10|l and JT7J). 
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FIG. 4: As Fig. [3] for the power-law potential. 



Equation (f2T|) indicates that, for a given C, our models 
approximate a skating solution at sufficently early times 
(small Q^). We can demonstrate this explicitly by solving 
equation ([5]) for the case A = 0, giving 

(l-w)(l+w) \ =k, (22) 
(1 — 11^) 

where k is a constant. At first glance, it is not obvi- 
ous how equations ([22)) and (|T7J|) could correspond to the 
same behavior for w(fl^). However, if we take the limit 
of small 0,$ in equation (|10|) . as indicated by equation 
(|2"Tj) . and take 1 + w <C 1 in equation ([2"12j) (since it is 
only in this limit that equation (fT0|) will be valid), we see 
that the two expressions are. in fact, equivalent, with k 
in equation ([22} given by k = (2/3)A§C 2 . 

We explore this equivalence in Fig. [5j Here we plot 
several w(a) curves for the linear potential, along with 
the evolution of w for the corresponding pure skating po- 
tential. As expected, the two sets of trajectories coincide 
at early times but diverge at late times. However, since 
it is the latter which is the epoch of interest observation- 
ally, the approximation developed here cannot usefully 
be replaced by the solution appropriate to a constant 
potential. 



IV. DISCUSSION 

We have extended earlier work on nearly flat poten- 
tials to include the cases where <j> 7^ initially. In these 
cases, the single trajectory for (1 + w)/Aq as a function of 
Q^, derived in Ref. [25j . becomes a family of trajectories 
parametrized by a constant C that depends on the initial 
value of <f>. All such trajectories show roughly the same 
behavior: an initial period of freezing behavior, for which 
w decreases with time, followed by a thawing period, for 
which w increases. However, this thawing period can 
take place in the future, so these models allow for purely 
freezing behavior up to the present. Note that behav- 
ior of this type (with freezing followed by thawing) was 
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FIG. 5: The equation of state parameter, w, as a function of 
the scale factor, o, for the linear potential (equation [18} with 
a — 0.1 and the indicated values of C (solid curves), along 
with the corresponding models with a = (dashed curves). 



also noted in the Monte Carlo simulations of Hutcrer and 



Peiris [44j , although their space of sampled models differs 
significantly from the class of models explored here. 

Our results differ from the earlier work in Refs. (32[ 
and [36j . who also examined freezing models. However, 
these papers assumed a model that began on a tracking 
trajectory, with w initially constant and far from —1, and 
then examined the freezing evolution as the scalar field 
began to influence the expansion. The models discussed 
here assume that the potential is sufficiently flat that w 
is already close to —1 at late times, so we expect very 
different evolution between these two cases. 

The results presented here can be easily generalized, as 
in Ref. (27j . to phantom models with a negative kinetic 
term. The result is simply equation (|10[) multiplied by 
— 1 on the right-hand side. 

Obviously, our results apply only to a special set of 
quintessence potentials. However, they do provide an 
interesting set of restricted evolutionary paths for such 
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potentials, and potentials satisfying the slow-roll condi- 
tions are a "natural" way to produce the observed value 
of w close to —1 today. One might argue, in that regard, 
that the freezing models discussed here are less natural 
than the thawing models discussed in Ref. [25j , This 
is a valid criticism, since in the former case one must 
tune the initial value of 4> to give w sufficiently close 
to —1 at present. Furthermore, the evolution discussed 
here cannot be extrapolated arbitrarily far into the past, 
since this would result in an early universe dominated by 
scalar held kinetic energy, in contradiction to observa- 
tions [45|. Thus, the slow roll conditions on V((f) cannot 
be satisfied at arbitrarily early times; the field must have 
evolved from a region with a different form for V(4>), giv- 



ing a different evolution for w(a) in the early universe. 
This means that freezing slow-roll models are necessarily 
more complicated than the thawing models considered 
in Ref. [25(, which do not have similar problems with 
their evolution at early times. However, in the absence 
of any compelling a priori models for the scalar field, it 
is worthwhile to consider all reasonable possibilities. 
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